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Abstract 

i— i A rotating superfluid such as a Bose-Einstein condensate is usually described by the Gross- 

i-Eh Pitaevskii (GP) model. An important issue is to determine from this model the properties of the 

quantized vortices that a superfluid nucleates when set into rotation. In this paper we address the 

S minimization of a two dimensional GP energy functional describing a rotating annular Bose-Einstein 
condensate. In a certain limit it is physically relevant to restrict the minimization to the Lowest- 
+Jj Landau-Level, that is the first eigenspace of the Ginzburg-Landau operator. Taking the particular 



^} MSC: 35Q55,47J30,76M23. PACS: 03.75.Hh, 47.32.-y, 47.37.+q. 

Keywords: Superfluidity, Gross-Pitaevskii energy, Bose-Einstein Condensates, Vortices. 

Contents 

(N 

> 1 Introduction [2] 

1.1 The GP energy and the LLL reduction [3] 

1.2 Main Analytical Results 

t~H 1.3 Discussion [8] 

2 Proofs QJ] 

2.1 The Minimization Problem QT] 

2.2 Infinite Number of Zeros [12] 

2.3 Construction of Critical Points [14] 

3 Numerical Simulations 1161 

3.1 Method Q1 

3.2 Results OH 



1 



2 



1 Introduction 

A Bose-Einstein condensate (BEC) is an object exhibiting quantum properties on a macroscopic scale. 
The effects of such phenomena as superfluidity can be observed experimentally on these systems, in 
particular the nucleation of quantized vortices in a rotating BEC (see for example the review |Fe2j 
or |MAH[ IMCWl IARV[ IRAV) IHCE|. IECH] for reports on original experiments). The most common 
description for a BEC is through a macroscopic wave function i/j : R d H> C which modulus squared will 
give the matter density profile. Here d = 2, 3 and in this paper we will focus on a situation where it is 
physically justified to consider a two-dimensional BEC. For a condensate at equilibrium in the rotating 
frame, this wave function should minimize the following Gross-Pitaevskii functional 

£ GP m = ! (\ |W> - in x ^\ 2 + (v(x) - ^\x\A + ~\<pA dx, (i.i) 



under the mass constraint 



H 2 = l. (1.2) 



Here x = (xi,X2) S R 2 , x 1 - = (—X2,xi), is the angular velocity at which the condensate is rotated 
around the axis perpendicular to the plane x\,X2 and G is a dimensionless coefficient characterizing the 
strength of atomic interactions. An essential feature of the experimental set-ups is the trapping potential, 
denoted by V in ( |1.1[ ). It is this potential that confines the atoms in a bounded region and prevents them 
from flying apart under the action of the centrifugal forces. For most experiments it takes the form (in 



the appropriate units that we have implicitly chosen in (1.1)) 

V{x) = \\x\\ (1.3) 
Such a potential sets a limit to the rotational speed 17 that one can impose to the condensate. Indeed, one 



can see that with V chosen as above the second term in ( 1.1 ) is not bounded below if ft > 1. Consequently 
the minimization problem makes no sense. Physically this corresponds to the fact that the centrifugal 
force overcomes the trapping force and drives the condensate out of the trap. 

To prevent such a singular behavior in the limit £1 — >• 1, Fetter [Fc2 proposed to use instead a potential 
of the form 

V(x)= 1 -\x\ 2 + ^\x\ i (1.4) 

which, at least theoretically, allows for arbitrarily large rotation speeds f2. Such a potential is a good 
approximation for those that have been used in the experiments ( |BSS[ ISBCj ) and this is the one we 
shall consider in this paper. 

An important experimental test for the superfluidity of a BEC is the observation of quantized 
vortices, which are described in Gross-Pitaevskii theory as zeros of the wave function carrying a positive 
topological degree or winding number (corresponding to the quantized phase circulation around the 
vortex). When a BEC is rotated, it will contain more and more vortices as the rotation speed is 
increased, until they become densely packed in the condensate when Q approaches 1. Vortices repel each 
other, and to minimize this effect they will arrange on regular hexagonal lattices when Q — > 1, for both 
the harmonic trap (1.3) and the 'anharmonic' one ( |1.4[ ). These Abrikosov lattices (named by analogy 



with the physics of type II superconductors, | Abr| ) have been observed experimentally, see e.g. |ARVj 



for experiments with the trap (1.3) and [BSSj for experiments with the trap (1.4). It should be noted 



that all vortices in these lattices are simply quantized, i.e. they have degree 1. 



A striking new feature of the trap (1.4) is the possibility to create annular condensates by increasing 



the rotation speed beyond the barrier that would be imposed by a weaker confinement. Indeed, for 
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sufficiently large rotation speeds, the centrifugal forces will dip a hole in the center of the condensate. 
Interesting questions then arise about the vortex structure of the condensate. For example : does 
the hexagonal vortex lattice survive the formation of the hole ? When the hole has just formed the 
answer is yes, see |FJS[ FB, IKBj for references in the physics literature and |CY) for a mathematically 
rigorous treatment. At very large rotation speeds however, the vortices disappear from the bulk of the 
condensate, resulting in a so-called giant vortex phase, a situation that we do not address here and for 
which we refer to |FJS[ IFBl 1KB! IKF] for physical references and |CRY[ iRou] for a mathematical analysis. 
The locations of vortices in the central hole of low matter density can also be investigated. It is 
commonly known (see e.g. fFJS| 1KB) IKTUj ) that there should be a phase circulation around this hole, 
indicating that vortices lie within it. In regions of high matter density, it is most favorable for the 
vortices to have degree 1, but it is far from obvious that it should also be the case in regions of low 
matter density. In particular, a question arising naturally |ABD[ IFJSj is wether the vortices in the 
central hole of an annular Bose-Einstcin condensate gather in a single multiply quantized vortex or form 
a pattern of simply quantized ones. 



In this paper we address the case of a two-dimensional annular Bose-Einstein condensate rotating 
in a potential of the form (1.4 1. We will study the formation of the annular condensate by taking the 
limit k — > of the problem, with the rotation speed properly scaled so as to capture the regime where 
the central hole is created under the influence of the centrifugal force. In such a parameter regime it 
is physically justified to consider a simplified variational problem that we describe in the next section. 
The very particular structure of the variational space considered allows to obtain interesting theoretical 
results. It also suggests a very natural numerical minimization method that we have used to confirm and 
extend our analytical results. The results we are going to present first appeared in a note intended for 
physicists |BR) to which we refer for further physical discussion (see also |AM| ). 



1.1 The GP energy and the LLL reduction 

We consider the GP energy functional 

£ GP m = J v Q |W - flte^l 8 + M a + ^l 4 ) H 2 + % M 4 ) dx, (1.5) 

to be minimized under the mass constraint 

M a = l. (1.6) 

The main idea of our analysis is to restrict the minimization of £ to the first eigenspace of the Ginzburg- 
Landau operator — (V — iClx^j , corresponding to the eigenvalue f2. This is the Lowest Landau Level 



(LLL), introduced in the context of Bose-Einstein condensation by Ho |Hoj . 

The use of such a simplification has originally been introduced for condensates rotating in traps of the 



form (1.3), i.e. for the functional 



£ GP [ip] = f Q |W - i£lxH\ 2 + l -^\x?W + % H 4 ) dx. (1.7) 

Restricting the minimization to the LLL is then justified in the limit Q — > 1. Indeed, the first term in the 
energy is at least equal to fl ~ 1 since f2 is the first eigenvalue of — (V — iflx- 1 ) . On the other hand, 
the second term will obviously be small in the same limit. Thus, minimizing the sum of the second and 
the third term one also obtains something much smaller than fi. 

This indicates that the energy must stay close to f2 when f2 — > 1. But the spectral gap between the first 



4 



and the second eigenvalue of the Ginzburg-Landau oper ator is equal to J7. To obtain an energy of order 
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fi the projection of a minimizer of the functional (1.7) on the excited energy levels of — (V — itlx- 1 ) 
must then be negligible. 

This formal justification has been the basis of several works on harmonically trapped BECs in the regime 
fl — > 1, e.g. jABD, CKR. WBP . A rigorous justification has been given in [AB2J. For other mathemat- 



ical developments on the functional (1.7) in the regime 51 —> 1 (Landau regime), we refer to |AB1[ RBN2 . 



The LLL consists (see e.g |LP| ) of functions of the form ip(z) = f(z)e n ' z ' / 2 where z is the complex 
variable xx + ixi and / is a holomorphic function. For an LLL function ip, the Gross-Pitaevskii energy 



(1.5) reduces to (we denote dz = dx\dx2) 



r ;i 'H = o : / (('^N 2 + |n 4 ) + '),/:. 

In this paper we will rather consider functions of the form i()(z) = /(z)e~' z ' I 2 (thus we effectively restrict 
the minimization to the first eigenspace of — (V — ix^j ) to have a functional space independent of ft. 
This approximation is harmless because, as we will see, our analysis is concerned with rotation speeds 
close to 1. As for the energy, the difference is very small as showed by |AB2( Lemma 3.1] : with if> of the 
form ip( z ) — /(^) e_ ' z ' ^ 2 with / holomorphic we have 

£ gp [tJj] = n + £ LLL [4>} (1.8) 

£ LLL [ip] = ^^(l-fi)|z| 2 + ^| 4 ^| 2 + ^|V,| 4 )rfz, (1.9) 

so that we will minimize the energy £ LLL amongst functions of the form ip(z) = f(z)e~\ z \ I" 1 with / 
holomorphic, under the mass constraint J \ip\ 2 = 1. For this purpose we will use the mathematical 



framework introduced in |ABN11 |A"B N2 for the study of the functional (1.7). See also |Nie2] for the 
corresponding dynamical model. 

Note that the validity of such a reduction could be investigated with the tools of |AB2) (as well as the 
validity of the reduction to a two-dimensional model), but in this paper we shall be c oncerned only 



with the study of the reduced energy functional £ LLL . We will provide in Subsection [L3j an a posteriori 
criterion for the validity of the reduction. 

We now describe the parameter regime that we shall consider. Let us define a small parameter 

£ = fc 1 / 3 (1.10) 

corresponding to a small anharmonicity regime and study the asymptotics of the problem as e — > 0. We 
take satisfying 

n = i + /3k 2/3 (i.ii) 

and will consider j3 and G as fixed. This choice will lead to a functional with bounded coefficients (see 



(1.20) below). In this regime the radius of the condensate is not bounded as a function of e. We thus 



rescale distances by making the change of variables 

^{z)=s- 1/2 ^{e- 1/2 z) (1.12) 

and for every <p we define the function 

f(z) = 0(z)el*l 2 / 2£ . (1.13) 
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By definition of the LLL, / belongs to the Fock-Bargmann space [Bar J 

5 e := j/ holomorphic ,J \ f\ 2 e^ 2/£ dz < ocj . 
The space 5" 6 is a Hilbert space for the scalar product 

(f,g) := / /(%(z) e - N2 /^z. 



(1.14) 



(1.15) 



The point of introducing such a space is that the orthogonal projector from L 2 ^C,e l z ' ^dz^j onto 3? e 
is explicitly known [Marl IFolj and called the Szego projector: 



n £ (<?)(*) 



us 



z'/e-\z'\ 2 /e 



g{z)dz. 



We will also use the spaces 



S| := j/ holomorphic ,J( 1 + NT |/| 2 e~ |z|2/e dz < ocj 



that we equip with the norms 



(l + |z| 2 ) s |/| 2 e -W 2 /^z 



1/2 



After the change of scale (1.12) and the change of function (1.13), the energy becomes 

£ LLL [iP] =eT LLL [f] 

with 



-/3|z| 2 + i|z| 4 )|/| 2 e^ + ||/| 4 e- 2 N 2 /^ dz 



and we thus have to minimize (1.20) under the mass constraint 

11/11*. = 1- 



(1.16) 
(1.17) 

(1.18) 

(1.19) 
(1.20) 

(1.21) 



It is this minimization problem that we will analyze in the sequel, namely we look at the problem 

F LLL := inf { J" LLL [/] \f€dl, H/ll*. = 1} • (1-22) 

In particular we will see that this problem admits a solution, i.e. the infimum above is actually a 
minimum. We are able, using the Szego projector (1.16) to derive the Euler-Lagrange equation satisfied 
by any solution of the minimization problem. Using this equation and the very particular structure of the 
Fock-Bargmann space (in particular the constraint that / has to be analytic) we are able to derive that 
any minimizer must have an infinite number of zeros (vortices) if e is small enough. We also construct 
almost critical points, functions that are solutions to the Euler-Lagrange equation up to a small remainder 
term. Evaluating the energy of these almost critical points yields an upper bound to the energy that we 
believe is optimal although we are not able to prove it. This upper bound will be useful to discuss the 
validity of the LLL reduction in Subsection |1.3| below. 
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1.2 Main Analytical Results 

Let us start by showing that the problem we want to consider is actually well-posed : 
Theorem 1.1 (Well-posedness and Euler-Lagrange equation). 

For any e > fixed, the problem (1.22) admits a solution in Any minimizer is a solution to 



-(3M E f + - (M 2 J + eM £ f) + GIl £ (e~^ 2 / £ \f\ 2 f) = (if 



(1.23) 



where [i is the Lagrange multiplier coming from the mass constraint and M e is the operator defined by 



MJ = ed z (zf) 



Alternatively, the Euler-Lagrange equation (1.23) may be written 

(-13 + \e)Il E (\zff) + \li e (|z| 2 n e (k| 2 /)) + Gn £ (- |z|2/£ |/| 2 /) = M/ 



or 



-ISMef + ) {M 2 £ f + sMj) + %f{ed z )[f 2 (z/2)] = [if, 



(1.24) 

(1.25) 
(1.26) 



where the operator f(ed z ) is defined as 



f(ed z 



+ 00 

k=0 



if f( z ) = T,t=o a kZ k 



Next we show that a minimizer of ( 1.22 ) cannot have a finite number of zeros if e is small enough. In 
particular, the minimum cannot be achieved by a polynomial. 

Theorem 1.2 (Infinite number of zeros). 

If e is sufficiently small, any solution to the minimization problem \1.22 ) has an infinite number of zeros. 



We will discuss this result in further details in Subsection 1.3 and now describe the procedure we 
follow to construct approximate solutions to the Euler-Lagrange equation (1.25). 

As is the case for harmonically trapped condensates in fast rotation, we expect that in the range of 
parameters we explore the scales of the problem will decouple. Namely we expect any minimizer of ( 1.20 1 
to be of the form cm where u varies on the scale of the vortex pattern (which is small compared to the 
size of the condensate) and a is a slow varying profile giving the general shape of the condensate. The 
kind of profile a we have in mind is a Thomas-Fermi distribution. More precisely, a should look like a 
minimizer of the energy ( |1.20 1 without the holomorphy constraint. Such a function has compact support 
and can thus not be in the LLL, but one can approach it by an LLL function which is an almost critical 
point for the energy (1.20). We proceed as follow. Let us introduce 

- ]z]2/2£ fr{z), fr{z) 



u T {z) 

where t = 17? + irj is any complex number and 

+00 



7 fc 



n e 27rr(n+l/2) 2 e (2n+l)7rii; 



(1.27) 



(1.28) 
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is the Jacobi 9 function (more precisely the 0o function according to the classification of |Niel| ). It has 
the property (see |Chaj for more details) 

e(v + k + It, t) = (-l) k+l e- 2iwlv e- MT Q(v, r). (1.29) 

The interest of introducing such functions is twofold. Firstly it is known |Cha| that any function v whose 
modulus is periodic over the lattice yf^Z © ^^Zt, vanishes exactly on the points of the lattice with 

simple zeros and such that g — ue' z ' I 2e is holomorphic must be proportional to u T . Secondly, the function 
f T is a solution to the Abrikosov problem (see e.g. |Abrj and [ABN21 Section 4] for a detailed discussion) 

n £ (|M 2 e-| Z ! 2 / £ / T ) = X T f T , With \ T =j \u T \ 2 b(T), (1.30) 

and 

6(r ) = J±lt_ = e -'l^-'l a /n. (L31) 

GfKI 2 ) Mez 



Here we denote j-v the average of a periodic function v. Equation (1.301 is similar to (1.251 without 
the potential term and with p = A T so that one can expect to obtain a solution of (1.25) by a slight 
modification of f T . We refer to |ABN1[ IABN2j and references therein for details on the quantity b(r). 
Let us just mention that it is minimum {b{r) ~ 1.16) for r = e 2l7r / 3 , which corresponds to a hexagonal 
lattice. From now on we take 

r = e 2 "/ 3 

and denote b = b(e 2l7r/3 ). 

To obtain our approximate solutions we first multiply f T by a slow varying profile a with compact support 
and J R2 \a\ 2 = 1. This trial state would correspond to having a regular hexagonal lattice of vortices in the 
bulk of the condensate (that is, the support of a). To obtain an admissible test function for the problem 
( 1.22[ ) we project af T onto the Fock-Bargmann space and normalize it : 



Accordingly 

u a T = f ayT e ' I / . (1.33) 

As for the choice of the appropriate a, it turns out to be related to the following Thomas- Fermi energy 
functional : 

^[P] = I ((-/3M 2 + JN 4 ) P+ ^P 2 ) dz. (1.34) 



The associated minimization problem is 



E^ F := inf \sr\p] , 9 > 0, / \z\ 2 p < +00, [ p=l 



(1.35) 



Its solution is unique and given as 



Pb ■= max f ^ ^ ii-L, Oj (1.36) 

where the chemical potential p[ F appears as a Lagrange multiplier in the equation associated to the 



minimization (1.35) ans is thus fixed by the constraint J K2 pJ F = 1. Note that there exists a critical 
value (3 C such that the support of such a function becomes an annulus when /3 > /3 C (see the discussion 
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in Section 1.3 1. This is in contrast with the situation where k — (purely quadratic trap) where the 



support of the TF minimizer is always a disc. 

We have the following result : 

Theorem 1.3 (Almost critical points). 

Let a be a profile satisfying 



\a(z)\ 2 = P r (1.37) 



and f aT be the associated function via formula (1.32). There holds 



(-/3 + i e )n £ (|z| 2 / Q ^) + Jn £ (|z| 2 n e (|z| 2 / Q ,r)) + Gn e ( e -i*i 3 / £ |/ a , T | 2 / Q! r) = ^ F U,r + Re (1-38) 

where \\R e \\$ e < Ce 1 / 4 . Moreover 

J" LLL [/ Q , T ] = jf ((-/3|z| 2 + \\z\^ |a| 2 + ^|a| 4 ) dz + 0(e^) (1.39) 



and thus 



F LLL < J7LLL[yr ] = ^TF + ( £ l/4)_ ^ 4Q ) 



1.3 Discussion 

In this subsection we discuss some of the physical insights that one can obtain from our theorems, 
and some of the questions that can not be answered analytically, for which we will rely on numerical 
simulations. 



Let us first comment a little bit more the ideas behind the construction leading to Theorem |1.3| It 
is known from the experiments that in the regime we consider, the condensate should contain a very 



large number of vortices. Theorem 1.2 also suggests to use a trial function containing as many vortices 
as possible, although it does not give information on their locations. On the other hand the trapping 
potential will force the condensate to live essentially in a bounded region, the complement of which will 
be a region of very low matter density. It is natural to think that, at least in the bulk, the vortices will 
form a regular hexagonal lattice. This is widely observed in experiments and numerical simulations. 
Combining these two requirements, we arrive at a trial function of the form caf T where f T is defined in 



(1.27), a is a function with compact support and c a normalization factor. Such a function is of course 
not in the Fock-Bargmann space : it has compact support and thus cannot be holomorphic. To obtain 
an admissible trial state we project it onto $ e using the Szego projector. Computing the energy of such 



a trial function consists of essentially two steps (see Subsection 2.3 for the detailed proof). First, using 



a remarkable result from IABN2] we have (see Lemma 2.1 below for a precise statement) 



n e (a/r) ~ ufr 

when e — > 0. Thus the projection onto the Fock-Bargmann space has no effect on the energy to leading 
order and we obtain 



T 



LLL 



y a , T ] ~ \mafr)\\t J c (-/%l 2 + |M 4 ) M 2 IM 2 e- |z|2/e ^ 

+ l|n £ (a/ T )|| Se 4 ^|a| 4 |/ r | 4 e- 2 N 2 /^z. 
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Now, |/ T |e l z l 2 / 2£ = |u r | is periodic over a lattice whose period is very small (oc e 1 / 2 ), whereas a is chosen 
with a support of fixed size. Therefore one should expect a homogenization effect leading to 



^ LLL [f a , T ] ~ ||n E (a/ T )||- 2 ^ K| 2 J (-f3\z\ 2 + \\z\^ \a\ 2 dz + \\n E (af T )\\^j \u T \^ J 



\a\ 4 dz. 



Computing the norm of H e (af T ) uses the same ideas and leads to 

1/2 



|n e (a/ r )b e 



/l«r| 



Finally, using (1.31), 



z\*\\«? 



bG, 



Id 4 dz 



and one immediately sees that the optimal choice for a is the one we made in ( 1.37 ). One should note that 



the only contribution of the vortex structure to the leading order of the energy is through the coefficient 
b, which depends on the type of lattice we choose and is minimum for the hexagonal lattice. 



We do believe that this construction is optimal, but a proof of a lower bound matching (1.40) would 



probably require to show that the vortices of a true minimizer of the energy are located close to the 
sites of a regular lattice. This is a difficult question, linked to the crystallization problem, and remains a 
challenging open problem in more than one context, see [XBN2, SSj and references therein. 
We have however the easy but non-optimal lower bound 



e7 f <f 



LLL 



(1.41) 



where E{ is defined by taking b = 1 in (1.35). This follows by minimizing T [/] with respect to 



j e -M /2e j dropping the constraint that this function should be in the LLL. As b is a constant, this lower 
bound confirms that Ej F gives at least the order of magnitude of the energy. The gap between the 
upper and the lower bound lies in the coefficient b that takes into account the density modification due 
to the presence of a large number of vortices. Filling this gap would be a first step towards a rigorous 
proof of the onset of the Abrikosov lattice in this regime. 



Although we have no rigorous proof that the trial state (1.32) accurately describes the true state of 



affairs, it is useful for a physical discussion to give some details on this function. The analysis of a profile 



such as (1.36) has already been carried out in [FJSUABD] . so we only adapt and summarize their results. 
The critical rotation speed for the condensate to develop a central hole is 



tt, = 1 



3fc 2 6G 
8tt 



1/3 



(1.42) 



For subcritical velocities, the behavior of the condensate is not qualitatively different from that of a 
harmonically trapped condensate, so we focus on velocities fl > Q c . In our scaling (see ( |1.11[ )) this is 
equivalent to the requirement /? > (3&G) 1 / 3 (87r)~ 1 / 3 . Then the inner and outer radii of the condensate 
R± are given by the relations 



1? 



R 2 _ = 4/3, R\-R 2 _ = (24bG) 



1/3 



the chemical potential is 



tf 

Mb 



/36G\ 



2/3 



(1.43) 



(1.44) 
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and the energy is 

«--!(£f-'- 

In the case of a disc-shaped condensate, the order of magnitudes of the corresponding quantities is the 
same as above. 

This information allows to check a posteriori the validity of our analysis. First, the LLL reduction should 
be valid if the energy F LhL is much smaller than the spectral gap between the Lowe st La ndau Level and 



the first excited level. In our scaling, this gap is oc e 1 = k 1 / 3 , and thus from (1.40) and (1.45) we 
deduce that the reduction is valid if 

G< fc~ 1/2 , /3< fc- 1 * 76 . (1.46) 

One could probably show rigorously that the LLL reduction is valid under these conditions, working in 
the spirit of [AB2J. Note that we have considered (3 and G as fixed in this paper so a rigorous proof 
would require to track the dependence on this parameters of our remainder terms. At least one condition 
will appear for sure if one wants to use our upper bound : we have used a (elementary) homogenization 
argument which requires the number of vortices in the support of p F F to be large compared to 1. Given 



that the vortices lie on a lattice of period oc e 1 / 2 = fc 1 / 6 and using (1.43) we arrive at the condition 



fc<G. (1.47) 



It is interesting to remark (tracking the dependence on ft and G in the proof of Theorem 1.2) that 



under the above condition the number of vortices for the true minimizer should be infinite. We also note 
that, as we require both k <C G and G <C necessarily k -C 1, which justifies our study of a small 
anharmonicity regime. 

Finally, interesting questions arise about the vortex pattern displayed by a minimizer of the energy. 



If one believes that a minimizer is close to our trial state (1.32), then the property 

n e (a/ r ) ~ af T 



that holds true in L°° (amongst other topologies, see Lemma 2.1 below) indicates that the vortices 
should lie close to the sites of an hexagonal lattice in the region of significant matter density, that is the 
support of pJ F ■ Very little information however can be obtained about vortices lying in the region of 
low density. 

These 'invisible vortices' have nevertheless a contribution : as shown in [AB1|IA"BD| . they help create the 
average Thomas-Fermi density profile. Indeed, if the vortex pattern was regular in the whole complex 
plane, one would obtain a Gaussian density profile. More precisely it is to be expected that the regular 
vortex lattice will be distorted close to the edge of the condensate, resulting in invisible vortices that 
will help shape the condensate. An interesting question (see e.g. jABD) iFJSl IKTU] ) is then, when f2 is 
above f2 c and the condensate has an annular shape, what is the vortex distribution inside the hole ? In 
particular : do vortices gather in a central multi-quantized vortex ? 

A nice feature of the LLL regime is that the condensate is completely determined by the location of 
its vortices, and one can thus minimize the energy numerically by varying the locations of the zeros of 
the wave-function. This gives a direct access to the optimal vortex pattern and allows to spot precisely 
the 'invisible vortices'. A direct minimization of the Gross-Pitaevskii energy would not give such an 
information because one should a posteriori look for density dips in a region where the density is already 
very small. 

We have adapted the numerical method of [ABDJ to the present setting and focused on the case where 
the condensate is annular. We refer to Section [3] for details. In particular we find invisible vortices both 
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outside the outer edge of the condensate and in the central hole of low density, with multiply quantized 
vortices appearing at the center of the trap. We also find a good agreement between the optimal value 



of the energy numerically computed and the upper bound (1.40), thus giving another argument for the 
optimality of our construction. 

The rest of the paper is organized as follows : in Section [2] we prove our analytical results, the proof 
of each theorem occupying a subsection. In Section [3] we describe our numerical method and present our 
results. 



2 Proofs 

2.1 The Minimization Problem 

We prove Theorem |1.1| In this subsection, all parameters are considered as fixed. 



Proof of Theorem \l.l\ We first note (see ( 1.41 1) that the energy functional .F LLL is indeed bounded below 
in the minimization domain we have chosen (1.22). 

Our energy is defined only for / g $ 2 , so w e take the convention that J- LLL [f] = +oo if / <^ J- 2 . Let 



(/n)neN be a minimizing sequence for (1.22) and u n — f n e ' z ' l e . We note that 



MVl 2 < 



1/2 



M 4 Kf 



1/2 



N 4 KI : 



1/2 



because of the mass constraint. We deduce that 



^ LLL [/„] > \J c \z\ 4 \u n \ 2 dz - \0\ Qf \z\ A \u n \ 2 dz 



1/2 



and thus the sequence (|z| 2 w n )„ is bounded in L 2 whatever the sign of /?. This implies that (/«)« is 
bounded in 3f . In ABN2] it is proved by using the Bargmann transform that 3* is compactly embedded 
in Similarly g| is compactly embedded [j] in 3"s- We thus have (possibly after extraction of a 
subsequence) the convergence of /„ to some limit analytic function /, strongly in $1, strongly in 3e and 
weakly in $ 2 . We deduce that 



inf 



^M 4 ) l/«| 2 e- |z|2/£ > J c (-/3M 2 + JM 4 ) |/| 2 e^ |2/£ . 



and 



On the other hand, denoting 



l/ll 



\f(z)\ p e-\ z ^ £ dz 



i/p 



it is proved in |Car|, Theorem 4] (see also |ABN2[ Section 2.2]) that for any entire function / 

„l/2 

||/(e 'OIU? < ll/IU?, for any < p < q such that e * < 



1 The Bargmann transforms maps 5| onto the space (l — d^. + £' 2 ) 3 L 2 ( 



(2.1) 



(2.2) 
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where /(e *•) is the function mapping z to f(e t z). Taking t = ^f— , p = 2 and q = 4 we deduce that 
for any entire function / 

(£|/| 4 e- 2 N 2 /^y /4 <C||/|| ys . 

Thus our sequence (f n )n is also bounded in A\ and we can assume (after a possible further extraction) 
that the convergence of /„ to / holds also in the weak topology of A\. In particular, by convexity, there 
holds 



inf / \f n \ 4 e- 2 ^ 2 / e dz > [ \f\ 4 e- 2 ^ 2 / e dz 
n Jc Jc 



which concludes the proof that 



T LLL [f]<infT LLL [f n } 



and thus the proof that / minimizes under the mass constraint. 

We turn to the derivation of the Euler-Lagrange equation. The weak derivative of the functional F 1 ^ 1 " 
at / along g is given by 

^ LLL (/) -9 = J c ((-/?M 2 + \\z\ 4 ) />- |z|2/s + f|/| 3 /fle- aWV ') dz. 

Using an integration by parts and d s f = d s g — on the first term and the fact that ^l 6 (g) = 9 on the 
second term, we obtain (1.23). 



Equation (1.26) is obtained from (1.23) exactly as in |ABN2i Proposition 3.2] with some algebra on the 
non-linear term. To get (1.25) we use an integration by parts to show that M e f = n £ (|;z| 2 /). Then 
M e 2 / = n £ (|z| 2 n e (|^| 2 /)) and we get the result. 



□ 



2.2 Infinite Number of Zeros 



We now prove that any minimizer of (1.22) has an infinite number of zeros. The argument is by con- 
tradiction and in two steps, using the two equivalent forms of the Euler-Lagrange equation (1.23) and 



(1.26). 



Proof of Theorem \1.S\ 

Step 1. Suppose / has a finite number of zeros. Then one may write f(z) = P(z)e ip ^ where P is 
a polynomial and tp is a holomorphic function. Now / e J £ and the condition J c |/| 2 e - ' z ' l e dz < oo 
implies that Re(<p(z)) < \z\ 2 /(2e). It is well-known (see |Boaj for example) that a holomorphic function 
can satisfy this condition only if it is a polynomial of degree less than 2. Therefore we know that 



/(*) = P(z)e a ^ 



(2.3) 



and the integrability condition on / implies «2 < l/(2e)- Injecting (2.3) in (1.23) and comparing the 
exponential growth of the different terms of (1.23) as in ABN2 yields a\ = «2 = 0. So, if / has a finite 
number of zeros, it is a polynomial. 

Step 2. Now, suppose / is a polynomial of degree n 

n 

f(z) = Y i a k z h 

k=0 
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and inject this in (1.26). The term 



-0M e f + - (M £ 2 / + eMj) 



is a polynomial of degree n, by the very definition (1.24) of M £ . Therefore (1.26) implies that the term 

G 



^)[/ 2 (z/2)]=]>>^) fc [/ 2 (z/2)] 



fc=0 



is also of degree n. But (ed z ) [f 2 (z/2)] is of degree 2n — k, so that it must be that at = for any k < 



Then / is of the form 
with 



f(z) = a n z n , 



I 2 = Ue n+1 n\ 



because of the normalization of / in $ e . Injecting this a last time in (1.26) yields 

G (2n)! 



(n + l)e + 



2 - yu + \{2e 2 + ne 2 + n{n ~ l)e 2 ) = 0. 



2 7re2 2 "+ 1 (n!) 2 r ~ ' 4 
Using the improved Stirling formula |Robj we obtain a condition on n: 

Ge- 1 / 12 e 2 n 2 e 2 3ne 2 

(i + /3e> - fine + 



(2.4) 



2ney/ri 2 4 4 

We now need to bound the chemical potential /i : taking the J^-scalar product of each side of ( 1.23 ) with 
/ yields 



(2.5) 



Here we used the fact that the spectrum of —/3M e + | (M| + eM e ) is bounded below by — /3 2 + ^ — jg . 
This follows by noting that the spectrum of M e is constituted by the eigen value s (n + l)h with (non- 



normalized) eigenvectors z n . We use the upper bound on J-" LLL [/] of Theorem 1.3 (see Subsection 2.3 
the proof) to deduce 



for 



6 f 2.bG 
5 \~8^ 



2/3 o2 pe Is 2 a Ge-V 12 n 2 e 2 3ne 2 



(2.6) 



Minimizing the right-hand side of (2.6) with respect to n (taken as a continuous variable as it should be 
very large when e is small) for fixed e yields 



Ge 



-1/12 \ 2 / r ° 



2- 



-6/5 



(2.7) 



and we deduce 



G 2/3 > c ^4/5^2/5 _ ^ G 2/5 e -l/5 + ^2^ 

which is a contradiction if e is small enough. 



□ 
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2.3 Construction of Critical Points 



The proof of Theorem 1.3 requires the following result taken from |ABN2| . that we quote for convenience. 
Here we abuse notations by taking the norm of non-holomorphic functions. 

Lemma 2.1 ( |ABN2 ] Estimates for f a . T ). 

Let 7 G C 0:1 / 2 (C,C) have compact support in Br . If g G 3 £ for some s, then H E (-jg) G $ £ for any 
s' G K and 

l|n £ (7ff) -7511,.' < Csy^hllco^hh^ 4 - (2-8) 



Ifv = e~\ z \ 2 / £ g G LP(B Ro ) for 1 < p < +oo tfien 

|(n £ (75)-75)e- |z|2/e || ip <C p , ito || 7 ||co.v 2 |l«IU P (^ )- 

We now provide the 

Proof of Theorem \1.S\ Let a G C 0,1 / 2 (C,C) have a fixed compact support and satisfy 



lal a = l. 



We first claim that for any v e M 

{-P + ^)n £ (M 2 / Q ,r) + *n £ (|z| 2 n e (|z| 2 /^ T )) + Gn £ (e~^ 2 ^\f a , T \ 2 f a , T ) - vf a , T 



(2.9) 



(2.10) 



||n £ (a/ T )||^n £ -(3\z\ z + -|z| 4 + Gb\a\ 2 _ „ a / T + R e , (2.11) 



2.1 



with ||i? e ||{f E < Ce 1 / 4 . The proof consists in a repeated use of Lemma 
Let p be a polynomial. We estimate, using Cauchy-Schwarz and (2.8| applied with g = f T and 7 = a, 

(\z\ 2 af T ~\z\ 2 Tl £ (af T ))pe-M 2 ^dz 



< c\\phM e {af T ) - af T \\n < ce 1/4 lbll».- 



Thus, by density of polynomial in ^ e 

||n e (k| 2 n £ (a/ T )) - U £ (\z\ 2 af T )\\ dE < Ce 1 ^. 



(2.12) 



Similarly, applying (2.8) with 7 = a\z\ 2 

f (\z\ 4 af T - \z\ 2 Tl £ (a\z\ 2 f T ))pe-^ 2 ^dz 
Jc 



< C\\phM e {a\A 2 fr) ~ a\z\ 2 f T \\$l < Cs^WpW 



and 



(|z| 2 n e (a|z| 2 / r ) - \z\ 2 n £ (\z\ 2 n £ (af T )))pe^ 2 ^dz 



< \\ P \\ Se \\lL £ (a\z\ 2 f T -\z\ 2 Tl £ (af T ))\\^ 

< \\ P hJ\af T -n e (af T )\U<CeV A \\ph e 



which implies 



|n e (|z| 2 n e (|z| 2 n e (a/ r ))) -n £ (klVr)|| 5e < Ce l ' A . 



(2.13) 



15 



We turn to the non linear term in the equation. Using (1.30) we have for any polynomial p 



&/K| 2 / \a\ 2 aTl e (e-W 2 ^\f T \ 2 f T )pe-W 2 ^dz = J \a\ 2 af T pe-^ 2 ^ dz 

= [ \a\ 2 ae-^ 2 ^\f T \ 2 f T pe~^ 2 ^dz 
Jc 

+ 0{e 1/4 ). (2.14) 



The second line is a consequence of (2.8) applied with 7 = \a\ 2 a and p = g. On the other hand, using 
Lemma |2.1| again 



J (e- |z|2/e |n £ ( a / T )| 2 n £ (a/ T ) - \a\ 2 ae-^ 2 ^\f T \ 2 f T ) pe'^^dz 



< W S J|n £ (a/ T )|| 5 J||a| 2 |^| 2 -|n e (a/ T )| 2 e-l z l 2 / e | 



+ l|Pb.l|a|li-||«r||i-||a/r ^ n £ (a/ T )|b £ < e 1/4 ||p|k 



and thus 



Recall that 



We claim that 



n e (e-l z l 2 / £ |n £ ( a / r )| 2 n £ (a/ T ) - b (J |u T | 2 ) \a\ 2 af T 

n e (a/ T ) 



(2.15) 



fa,r 



\Tl e (af T )\\^ 
\Tl e (af r )\\l = -f\u T \ 2 (l + O^) 



(2.16) 



Indeed, from Lemma 12. II we have 



l|n e (a/ T )HL = / MV| 2 + o( e V4). 

On the other hand |u r | is periodic over a lattice of period proportional to e 1 / 2 , thus for any 7 € C°°(C) 



l7l 9 KI 9 -/KI 2 jfl7l s 



< C||7|| c o (C ) 



and 



|7| 2 K| 9 -/K| a / M 2 <C £ 1/2 ||7llc 1( C)- 
Using a G C* ' 1 / 2 and interpolating between the two previous estimates we obtain 



\a\ 2 \uJ 2 = l M 2 \a\ 2 



(2.17) 



and thus ( |2.16[ ), recalling ( |2.10| . Gatherin g (|2.12[ ),( |2.13| ,( |2 15| > and ( |2.16[ ) we deduce that \2.ll\ holds. 
There only remains to choose a satisfying (1.37) and v — Ji^ to have 



-P\z\ 2 + ^\z\ 4 + Gb\a\ 2 -,<)o = 



1G 



and thus (1.38). 



The proof of the energy result (1.39) uses the same tools. By definition 



^ LLL [/ Q , T ] = ||n £ ( a / r )ll 5e 2 £ (-/?NI 2 + ^I 4 ) |n e (a/ T )| 



\-\*\ 2 /e 



+ l|H«(«/r)||^ I |n £ (a/ T )| 4 e- 2|z| /e 



working as above, with the L p estimate (2.9) instead of (2.8) we deduce 



^[f a , T ] = \\M^fr)\\^ j(-^\z\ 2 + \\z\^ \au T \ 2 + \\n e {afr)\\^ j c \au T f + 0{e 1 ' A ). (2.18) 
Next, using that |it T | is periodic over a lattice of period proportional to e 1 / 2 as above (see the proof of 



(2.17)), we obtain 



^ LLL [/ Q , T ] = \\IiMfr)\\^j\Ur\ 2 J^P\z\ 2 + \\z\^\aY 

+ \\U e (af T )\\^J K| 4 ^ £ |a| 4 + Oie 1 ^). 



There only remains to recall (1.31) and (2.16) to conclude that (1.39) holds. 



(2.19) 



□ 



3 Numerical Simulations 
3.1 Method 

We want to find a numerical approximation of the minimizer 6 of £ LLL in the LLL. We write 

<f>(z) = P(z)e- n ^ 2 ^ 2 (3.1) 

with P a holomorphic function. Since polynomials are dense in it is reasonable to fix an integer n 
and to restrict the analysis to functions (j> where P is a polynomial of degree less than n (see [ABN2, 
Section 6] for mathematical results on the validity of this approximation). We write our trial functions 
as in [XBT3] : 



z)=A\{(z-z. ] )e-^/\ 

3 = 1 



(3.2) 



where 



A 



-1/2 



L 2 



is the normalization factor. We will numerically vary the locations Zj of the vortices. An alternative 
method (used for example in [CKR for a harmonically trapped condensate) would be to take 



-Q\z\ 2 /2 



and vary the coefficients bj . The interest of our approach is to give a direct access to the exact repartition 
of vortices, whereas the alternative method would require to compute the roots of a polynomial of degree 
n, which is a delicate task for large n. In particular, varying the coefficients could probably not give the 
precise locations of invisible vortices located in regions of low matter density. Note also that there exist 
some good methods to directly minimize the GP energy (see |DH[ IDK] for recent developments), but 
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they probably could not give access to the locations of invisible vortices neither. Indeed, one would have 
to spot zeros or phase singularities in a region where the density is very small. 

The point here is that these vortices, although invisible, have a crucial influence on both the energy and 
the density profile of the condensate. One can be skeptical about the existence of density dips in a region 
where there is so to say no matter at all. But a vortex is also carrying a phase circulation (or superfluid 
current) and thus one can interpret the presence of invisible vortices as the existence of currents in the 
condensate that are equivalent to those generated by vortices in the low density region. Taking invisible 
vortex into account with a trial function of the form (3.2) is thus a way to evaluate more precisely the 
effect of superfluid currents in the condensate. 

We used a conjugate gradient method with a Goldstein and Price line-search. All the functions whose 
integrals have to be computed are of the form Polynomial x Gaussian. We thus used the Gauss-Hermite 
method, and we took enough Gauss points for the computations to be exact. This results in quite 
expensive calculations, but we have been able to numerically construct condensates with up to ~ 120 
vortices. 



3.2 Results 




Figure 1: Vortex structure and atomic density for G — 3, k = 10 4 , j3 — 1 (f2 — 1 = 2.2 10 3 ). There are 
67 vortices in total, the central vortex is constituted of 11 single vortices. 

We show in Fig[l] and Fig[2] typical examples of configurations we numerically computed. The 
qualitative features of the vortex patterns and atomic densities confirm our theoretical results and are 
in good agreement with existing theoretical and numerical studies |FJS1 1KB) . Note however that the 
numerics become quite intricate for large number of vortices, which accounts for the relative lack of 
symmetry of Fig(2j 

As was expected, the condensate develops a central hole and visible vortices are regularly distributed 
and densely packed in the annular region of significant atomic density. The hexagonal lattice is not 
clearly observed however, which seems to be due to the competition between the annular geometry of 
the condensate and the repulsion between vortices. The former tends to force vortices to form concentric 
circles while the latter favors the hexagonal lattice. This effect is apparent in several other simulations 
of rotating annular BECs, see t AMj IFJSj for 2D simulations and |AD1 IDan| for 3D ones. To observe 
clearly the hexagonal lattice one would need the annulus to be thicker so that it could contain more 
vortices. This would make the numerics all the more challenging. 

Our computations also show a distortion of the vortex pattern near the external radius of the 
condensate, resulting in invisible vortices in the exterior region of low atomic density as is the case for 
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10 - 



g - 




Figure 2: Vortex structure and atomic density for G = 3, k = 10 5 , j3 — 1 (f2 — 1 = 4.6 10 4 ). There are 
119 vortices in total, the central vortex is constituted of 20 single vortices. 




Figure 3: Two example of vortex configurations for G = 3, k = 10 4 ,/3 = 1(0 — 1 = 2.210 3 ), respectively 
with n — 60 and n = 67 vortices. 

harmonically trapped condensates [ABDJ. Some vortices also lie in the central hole as theoretically 
predicted and we can get information on their precise locations: we observe a distortion of the regular 
lattice near the inner radius of the condensate, resulting in isolated singly-quantized vortices encircling 
a central multiply quantized vortex. We computed configurations for which this central vortex has up 
to 20 units of circulation while there is a total of 32 units of circulation in the entire hole and 83 visible 
vortices. The number of vortices (both visible and invisible) increases with increasing f2 or /3, but since 
our scaling does not allow to explore a large domain of fl when G is fixed, we mainly varied k. The total 
vorticity of the system increases with decreasing k. 

All vortices do not have the same contribution to the energy: as n increases, the vortex pattern 
in the annular region of significant atomic density remains the same up to possible rotations, with the 
additional vortices first gathering in the central vortex, then constituting the distorted lattice near the 
inner boundary of the condensate and finally occupying the distorted sites beyond the external radius. 
With increasing n (see Figj4]), the energy reaches a first plateau when the central figure is formed, 
constituted of the visible vortices and the vortices in the central hole. A second plateau is reached 
when enough distorted sites beyond the external radius are occupied. For example, for the parameters 
corresponding to Fig{l]and Fig|2|(fc = 10" 4 , G = 3, = 1), the energy varies by ~ ±10 6 in relative value 
when n is increased from 60 to 67 (Figp| and the atomic density does not vary significantly. We find a 
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goo d agreement between our numerical results and the upper bound to the energy computed in Theorem 
1.3 The energy of our analytical trial function typically differs by ~ 10 -3 to ~ 10~ 2 in relative value from 



the energy numerically computed. The agreement between the radii theoretically predicted (1.431 and 



extracted from the numerical simulations is typically a little worse but still satisfactory (relative errors 
ranging from 10~ 2 to 1CP 1 ). Again, obtaining a better agreement would probably require to consider a 
thicker annulus containing more vortices. This would result in more difficult simulations. 



Figure 4: Minimum energy as a function of the number of vortices in the trial wave function (G = 3, k 
HT 4 ,/3 = 1). 
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